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Abstract 

We describe some recent applications of Tsallis statistics in fully 
developed hydrodynamic turbulence and in high energy physics. For 
many of these applications nonextensive properties arise from spatial 
fluctuations of the temperature or the energy dissipation rate. The 
entropic index q is related to the relative magnitude of these fluctu- 
ations. We concentrate on a recently derived formula for the energy 
dependence of q that is experimentally verified by fits of cross sections 
in e~^e~ annihilation experiments. Evaluating this formula for much 
smaller energies E of the order of the recombination temperature, one 
obtains the correct order of magnitude of the fluctuations of the cos- 
mic microwave background. Evaluating it for — > oo, one obtains 
(under certain assumptions) possible evidence for the existence of 6 
compactified dimensions, as predicted by superstring theory. 



1 Introduction 



Nonextensive statistical mechanics as introduced by Tsallis in 1988 and 
further developed by many others is a useful generalization of ordinary 

statistical mechanics, with potential applications not only for equilibrium 
systems but also for nonequilibrium systems with a stationary state. It is 
based on the extremization of the Tsallis entropies 



(1) 



rather than the Shannon entropy Si, which is contained as a special case 
in this more general setting. The Pi are the probabilities of the physical 
microstates, and q is the entropic index. 

Given two independent subsystems I and II with probabilities and pf^, 
respectively, the Tsallis entropy of the composed system I+II (with proba- 
bilities Pij^^^ = pIpY) satisfies 

Sl-^'' = Sl + Sl^-{q-l)SlSl^. (2) 

Hence the system is nonextensive, i.e. there is additivity of the entropy for 
q = 1 only. 

Broadly speaking, the formalism with q ^ 1 has so far been observed to be 
relevant for three different classes of systems. The first class is systems with 
long-range interactions (e.g. the second one is multifractal systems 

(e.g. |T^-[0), and the third one is systems with fluctuations of temperature 
or energy dissipation (e.g. ||T6[|- [^ ). 



In this paper we describe some interesting new physical applications of the 
nonextensive approach. Our examples mainly seem to be lying in the 3rd 
class. One example is fully developed hydrodynamic turbulence |P!^-P5|, 
another one the statistical behaviour of particles produced in high energy 
e"'"e~ collider experiments [^-[^. We will also deal with possible aspects 
of Tsallis statistics in the early universe [^-|^, and in particular point out 



a connection with the fluctuations of the cosmic microwave background. In 
the last section we will show that (under certain assumptions on the orign of 
the temperature fluctuations) the e"'"e~ scattering data may provide exper- 
imental evidence for the existence of compactified dimensions, as predicted 
by superstring theory. 
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2 Fluctuations of temperature and Tsallis sta- 
tistics 



Let us here sketch the way in which fluctuations can lead to nonextensive 



behaviour. This is a relatively new approach to nonextensitivity |1T6|-||20||, 
which is in particular relevant for nonequilibrium systems with a station- 
ary state. The fluctuations may be either local temperature fluctuations or 
fluctuations of an effective friction coefficient (see ||3^, ^ for some related 
nonequilibrium systems) . 

Consider a system of ordinary statistical mechanics with Hamiltonian H. 
Tsallis statistics with g > 1 can arise from this ordinary Hamiltonian if the 
inverse temperature (3 is fluctuating. For example, we may think of these 
fluctuations as spatial fluctuations. From the integral representation of the 
gamma function one can easily derive the formula 0, ^ 

1 

(1 + (g - l)P,H)-— = / e-^^/(/3)ci/5, (3) 

Jo 



where 



is the probability density of the (or gamma) distribution. The above for- 
mula is valid for arbitrary Hamiltonians H und thus quite important. The 
left-hand side of eq. (^) is just the generalized Boltzmann factor emerging 
out of nonextensive statistical mechanics. It can directly be obtained by 
extremizing Sq. The right-hand side is a weighted average over Boltzmann 
factors of ordinary statistical mechanics. If we consider a nonequilibrim sys- 
tem with fluctuating j3 variables, then the generalized distribution functions 
of nonextensive statistical mechanics are a consequence of integrating over 
all possible fluctuating inverse temperatures /?, provided P is distributed. 



Wilk et al. [17| have also suggested an extension of this approach which 
is valid for g < 1. However, in that case the distribution function /(/?) is 
more complicated and different from the simple form (^. In particular, for 
g < 1 one has /(/?) = f{P,H), i.e. the probability density of the fluctuating 
temperature also depends on the Hamiltonian H. This case seems physically 
less straightforward, hence we restrict ourselves to g > 1 in the following. 
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For g > 1, the distribution obtained in eq. (H) is a universal distribution 



that occurs in many very common circumstances For example, it arises 
if 13 is the sum of squares of n independent Gaussian random variables Xj 
with average 

n 

p=j:x!, (5) 

i=i 

with 

n = (6) 
q-1 ^ ^ 

Hence one expects Tsallis statistics with g > 1 to be relevant in many appli- 
cations where n (approximately) independent and (approximately) Gaussian 
variables contribute to the variable (3. Of course, other distribution functions 
/(/?) can also be considered which may lead to other generalized statistics. 
Tsallis statistics is just one possibility, though a very relevant one. It is 
distinguished by generalized Khinchin axioms and formal invariance of most 
thermodynamic relations [0, ^. 

The constant Pq in eq. (H) is the average of the fluctuating P, 



oo 



(/?) := / PfimP = (7) 







The deviation of q from 1 can be related to the relative variance of (3. One 
can easily check that 

a i-M^ (8) 

3 Application to turbulent flows 



A lot of recent papers |18h|^| showed that extremizing Tsallis entropies is a 



very useful tool in fully developed turbulence, leading to predictions for scal- 
ing exponents and probability density functions that coincide quite precisely 
with experimental data. Different models were developed, some of them us- 
ing Tsallis statistics with a constant g < 1 within the multifractal model 
of turbulence others doing Tsallis statistics with a scale-dependent 



g > 1 directly in velocity space |^-|^. Apparently, the way Tsallis statis- 
tics can be used in turbulence is not unique, but often the predictions given 
are quite similar. 
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Here, let us follow the approach introduced in EM and experimentally 



verified in p5[. Consider turbulent velocity differences u := v{x + r) — v{x) 
at a certain scale r of the flow. Consider a nonlinear Langevin equation of 
the form 

ii = -^F{u) + aLit) (9) 

where F{u) = —-^V{u) is a nonlinear forcing and L{t) is Gaussian white 
noise. To be speciflc, let us assume that V{u) = C|up" is a power-law 
potential. The parameter f3 is deflned as 

P = ^. (10) 

For a = 1 we just obtain the ordinary Ornstein-Uhlenbeck process, for which 
P = •y/cr'^ is related to the temperature of the surrounding heat bath by 
P = 1/kT. Given some P and a, the stochastic differential equation (^ 
generates the stationary probability density 

p{u\P) = {CP)^exp {-PC\U\'-} . (11) 

So far this invariant density deviates from what is measured in turbulence 
experiments, so our model is yet too simple for turbulence. But now comes 
the important point. Let us assume that not only L{t) is a fluctuating quan- 
tity but also P = yjcy^. We assume that P changes spatially on a relatively 
large scale, i.e. the test particle flrst moves in a certain region with a given /3, 
then moves to another region with another P and so on. Within each region 
there is local equilibrium. If /? is a distributed random variable in the var- 
ious regions, with an n given by eq. (^, then for the marginal distributions 
•p{u) = J p{u\P)f{P)dp we obtain after a short calculation 

.si 1 
P{u) = 12) 

Zg{l + {q-l)pC\u\^'^) — 



where 



and 



2a , , 

q = 1 + ^— 14 

an + 1 

P = TT^ Po. (15) 



1 + 2a - g 
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The densities (|T^ extremize the Tsalhs entropies. Hence Tsalhs statistics 
is rigorously proved for this simple model with fluctuating friction forces 
pO| . Moreover, for a test particle moving in a turbulent flow the dynamics 
is physically plausible, since one knows that the dissipation fluctuates in 
turbulent flows. 

The densities p{u) as given by eq. (|1^) were just those that were recently 
successfully used as fits of experimentally measured turbulent densities p5| . 



They perfectly coincide with the experimental data (Fig. 1 in [^, Fig. 1 
in thus making this nonextensive model a very good model of turbu- 

lent statistics. The generalized Langevin equation yields a dynamical reason 
for Tsallis statistics. If one replaces the Gaussian white noise in the above 
stochastic differential equation by deterministic chaotic noise, one can also 
understand the (very small) asymmetry (skewness) of the observed distribu- 
tions ETl ra, 



4 Application to e+e annihilation 

Not only velocity distributions in hydrodynamics, but also momentum dis- 
tributions of particles produced in high-energy scattering experiments can 
be successfully described using nonextensive methods. The differential cross 
section measured in a collider experiment is proportional to the prob- 
ability density to observe particles with a certain momentum p. Bediaga 
et al. ||2^ looked at experimentally measured differential cross sections of 



c ■ u 



roo , I 2_ 

dx (l + {q-l)^x^ + u^ + m}) (16) 



transverse momenta pt in e'^e annihilation experiments and showed that 
the formula 

1 da 
a dpT 

yields an extremely good fit of the experimentally measured cross sections. 
Here x = pL/kT^, u = px/kTo and := mo/fcTo are the longitudinal mo- 
mentum, transverse momentum and mass of the particles in units of the 
Hagedorn temperature Tq, respectively, c is a constant related to the multi- 
plicity M, i.e. the number of charged particles produced. 

The Hagedorn temperature Tq is a kind of phase transition temperature 
of 'boiling' nuclear matter. This fundamental concept was introduced by 



Hagedorn in a seminal paper in 1965 ||3T|, |4^, These days the Hagedorn 



phase transition is seen in connection with the QCD phase transition from 
confined to non-confined states. 
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Bediaga et al. just wrote down eq. (|TB|) without deriving it. They ne- 
glected m0 and used three free fitting parameters for their fits, namely c, Tq, q. 
The question was if one can derive suitable equations for these fitting param- 
eters, thus proceeding from a mere fit to a theory. Although the complete 
theory taking into account all possible interactions within the hadronization 



cascade is out of reach, a first step was performed in |30|, where eq. (|T6|) was 
derived from first principles by g-generalizing Hagedorn's original approach. 
The tool for this are suitable factorizing nonextensive partition functions for 



fermions and bosons, actually of similar form as previously suggested in ||34 

One can then go a step further by considering suitable approximations 
and by making plausible assumptions on the energy dependences of various 
parameters, being guided by the existence of moments. The final result of 

e has 

M-p{u) (17) 



30| was that in good approximation one has 

I da 1 



a dpT To 



where 



p{u) = -^u^/^ {1 + {q - l)u)-^^^ (18) 



(19) 



To(g) = (i-|)pt (20) 

p1^ ^ 180 MeV (21) 

11 _ p-i?/i?0 

m = -TTT-WJ^ (22) 



9 + e-^/^o 
1 

2 



Eo ^ -mz = 45.6 GeV (23) 



/ p \5/ll 

M(E) « (^) . (24) 

Here E is the center of mass energy of the beam, and the function B{x, y) = 
T{x)T{y)/T{x + y) denotes the beta function. 

Formula (0) with p{u) given by eq. (plSj) , q{E) given by eq. (p2[) , To(g) 
given by eq. ( pUD and multiplicity M{E) given by eq. (|2^) turns out to 
very well reproduce the experimental results of cross sections for all energies 
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E in the experimentally accessible region 14... 161 GeV (see Fig. 3 in 1 3D 
which compares with experimental cross sections measured by the TASSO 
and DELPHI collaborations [|5|, ||). 

The theoretical progress is that eqs. (|17[)-(|^) provide concrete equations 
for the energy dependence of all parameters Tq, M,q used in the fits, and only 
two fundamental constants are left: One is p^, the energy scale of the QCD 
phase transition, and the other one Eq, the energy scale of the electroweak 
phase transition. 

The agreement of eq. {\n\} with the experimentally measured cross sections 
is remarkably good. Nonextensitivity seems to be the right approach to tackle 
this high energy problem, q is not constant but depends on the spatial scale 
r = h/E, just similar as it does in turbulence. Of course, an alternative for 
particle physicists is to do Monte-Carlo simulations. The simulated Monte- 
Carlo cross sections look very similar to our theoretical curves |^^. Thus 



the complexity inherent in the Monte- Carlo simulation effectively reduces to 
some model with nonextensive properties. 



5 Fluctuations of the cosmic microwave back- 
ground 

Many authors have already dealt with possible applications of nonextensive 
statistical mechanics in high energy physics. Work so far has concentrated 
on the solar neutrino problem |^-||5ll], quark gluon plasmas ||5^, heavy ion 
collisions [^, cosmic rays |54] 



microwave background [B^, p4 



e^e" annihilation 0, and the cosmic 
Useful bounds on Iq' — 1| in the early uni- 
verse were derived [^-[^. Here we want to deal with the cosmic microwave 
background from a new point of view. We will estimate the order of mag- 
nitude of the temperature fluctuations using our formula for q{E) from the 
previous section. 



Although formula (0) is only experimentally confirmed for energies in 
the region 14... 161 GeV, we can certainly also evaluate it for other energies 
as well. For E << Eq we may write exp{—E/EQ) ^ 1 — E/Eq, which yields 



q{E) 



1 E 
5^ 



{E « Eo). 



(25) 



In section 2 we saw that any g > 1 can be interpreted in terms of spatial 
fluctuations in (3 = 1/kT. It is now interesting to see what the typical size 
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of these fluctuations is if we evaluate our formulas at the temperature E ^ 1 
eV where the universe becomes transparent (the recombination temperature). 
One obtains from formulas (|8|), (1231), (|2^ 



/^^J^^2.10- (26) 



Here A(3 denotes the standard deviation of /3. Remarkably, we obtain a value 
that has the same order of magnitude as the experimentally measured value 
of the fluctuations of the cosmic microwave background : 



AT 

1-5 ■ 10"*^ < ^ < 10 ■ 10"^ (27) 

Of course, the correct order of magnitude may just be a random coincidence, 
but it may also not be! If true, then nonextensive statistical mechanics is 
probably relevant for the thermodynamic description of the early universe, 
due to fluctuations that were neglected so far. We expect that q is energy 
dependent and given by eq. (|2^ ) or something of similar order of magnitude. 
Our value of g — 1 ~ 4 ■ 10~^^ at the recombination temperature is certainly 
consistent with previous bounds derived, which were of the order |g — 1| < 
10-^ 



6 Compactified dimensions 

Another striking observation is related to the E ^ oo limit, where our for- 
mula (|2^) yields 

g(oo) = ^. (28) 

The largest possible energy E corresponds to the smallest possible scale 
r h/E. In a quantum gravity setting this is the Planck scale. It is now 
reasonable to assume that at the smallest possible scale r (described by the 
largest possible energy E) the number n of independent contributions Xj to 
the fluctuating /3 is given by the number d of spatial dimensions of the space. 
The reason is that the momenta whose kinetic energy define /3 can fluctuate 
in each space direction independently. At the smallest scales, there is (by 
definition of the smallest scale) no other structure and degree of freedom ex- 
cept for the various directions in (i- dimensional space. Using this (somewhat 
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speculative) argument we obtain from eq. (|^) 

2 2 2 

g — 1 = — = - = - (at smallest scales), (29) 
n a 9 

hence d = 9 at the smallest scales. 9 spatial dimensions mean that space-time 
is 10-dimensional. Remarkably, we get the number of space-time dimensions 
that are necessary for superstring theory to be formulated in a consistent 
way pB . 



The e+e" annihilation data actually already show a crossover from q ^ 1 
to g ~ 11/9 at the electroweak energy scale of about 100 GeV |2^, Should 
the saturation of q near the value 11/9 be confirmed by future experiments 
with larger center of mass energies E > 161 GeV, this could be interpreted 
as possible experimental evidence that the compactified extra dimensions of 
superstring theory start already becoming visible at the electroweak scale. 
This means, their diameter would be of the order of r ~ hc{100GeV)~^ , much 
larger than the Planck length. This is compatible with other theoretical 
arguments, which also suggest relatively large sizes of the extra dimensions 



||57|| . If this picture is correct, then the ultimate reason for the validity of 
Tsallis statistics in the early universe could be fluctuations originating from 
the Planck era, where q = 11/9. 
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